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Gelation in Aerosols; Non-Mean-Field Aggregation and Kinetics 

C.M. Sorensen and A. Chakrabarti 
Kansas State University 
Manhattan, Kansas 66506 

Abstract 

Nature has many examples of systems of particles suspended in a fluid phase; colloids when in a 
liquid, aerosols when in a gas. These systems are inherently unstable since if the particles can come 
together, van der Waals forces will keep them together. In this work we studied the aggregation kinetics 
of particulate systems, most often aerosols. The emphasis of our work was to study dense systems and 
systems that gel since previous work had not considered these. Our work obtained a number of significant 
discoveries and results. 

1. We found it useful to define three regimes of aggregation based on the relative distance between 
clusters: the cluster dilute regime where the nearest neighbor separation is very large compared to 
the cluster size, the intermediate regime where the nearest neighbor separation becomes 
comparable to the cluster size, and the cluster dense regime where the nearest neighbor separation 
is so small that the clusters travel with a nondiffusive, straight line motion between collisions. 
Previous work had concerned only the cluster dilute regime. 

2. We found that the mean field description of aggregation continues to work well as the system 
evolves into the intermediate and dense regimes. 

3. We found that aggregation in the cluster dense regime leads to the formation of superaggregates, 
i.e., aggregates with a fractal dimension of 2.6 over large length scales made up of smaller 
aggregates with a fractal dimension of 1.8 over their smaller length scales. Such superaggregates 
were found in the soot from combustion of heavily sooting fuels. 

4. The kinetics crosses over from the well known dilute case to the intermediate and dense cases as 
identified by an increasing aggregation kernel homogeneity. 

5. These results above allow us to describe gelation as the evolution of a particulate system from 
cluster dilute through the intermediate regime to the cluster dense where the system ultimately 
gels. 

6. Based on our understanding of gelation we, for the first time, caused an aerosol to gel. These 
aerosol gels have unique properties and a patent has been filed. 

7. Simulation of charged aerosol aggregation showed that the presence of a small fraction (10% or 
less) of charged monomers does not yield significant differences in the growth kinetics. 

8. We also studied how light scatters from these aggregates relative to scattering from simple 
spheres. We did this because our experiments use light scattering to probe the systems. We found 
heretofore undescribed patterns in the scattering as a function of angle which give insight into the 
nature of the scattering process. We also studied multiple scattering in dense particulate systems. 

Superaggregates in Dense Particulate Systems 

Flames 

We have shown that soot formed in laminar diffusion flames of heavily sooting fuels evolves through 
four distinct growth stages which give rise to four distinct aggregate fractal morphologies. Each 
successive stage grows upon the previous stage; hence aggregates of one morphology, parameterized by a 
fractal dimension, form larger superaggregates or gels of a different morphology with a different fractal 
dimension. These results were inferred from large and small angle static light scattering from the flames, 
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microphotography of the flames, and analysis of soot sampled from the flames. These results and the 
analysis were substantiated by comparison to computer simulations. The growth stages occur 
approximately over four successive orders of magnitude in aggregate size and involve either diffusion 
limited cluster aggregation or percolation in either three or two dimensions. Table 1 gives an overview 
our results, figures 1 to 4 show examples. 

TABLE L— THE REALM OF SUPERSOOT 


Stase 1 

Stase 2 
percolation or 

Stase 3 

Stase 4 

DLCA 

2d, DLCA 

2d percolation? 

restructuring? 

D^1.8 

-A D - 2.6 

-> D = 1.4 

-> D=1.9 

aggregates 

superagg(SA) 

supersuperagg(S S A) 

gel network 


~ 1 pm 

~ 20 pm 

~200 pm 




Figure 1. — Light scattering showing both the D = 2.6 and 1.8 morphologies, stage 2 and 1, of 
soot superaggregates in an acetylene/air diffusion flame. Also shown are soot aggregates 
extracted from the flame with these morphologies. 


NAS A/CR— 2008-2 1 5280 


2 



Figure 2. — Light scattering from a variety of flames some of which do not show superaggregates, some do. 
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Figure 3. — Example of 2-D stage 3 
soot. 


Figure 4. — Example of 
stage 4, continuous 
2-D percolation, soot. 


Simulations 

These results can be understood with the aid of our simulations of off-lattice cluster-cluster 
aggregation for diffusion-limited (DLCA), ballistic-limited, and reaction-limited cluster aggregation 
classes. We find that as the system evolves and becomes dense, the largest cluster develops a hybrid 
structure with mass fractal dimension Df= 2.6 over large length scales, while at smaller length scales, the 
early time dilute-limit fractal structure with D = 1.8 is frozen in, figure 5. Consistent with our experiments 
the largest cluster is thus an aggregate of smaller aggregates with a different fractal dimension, and we call 
it a “superaggregate.” The crossover length separating the two morphologies, which we call the critical 
radius of gyration, can be calculated based on a simple theory that assumes a monodisperse cluster size 
distribution. This agrees well with simulation results for DLCA. 

Chamber Explosions 

We use 4 and 17 liter chambers (fig. 6) to create dense, nanoparticle aerosols which will gel. The 
dense aerosol is created by rapid (explosive) reaction of gases such as C 2 H 2 or SiH 4 with an oxidizer (e.g., 
0 2 ). From digitized pictures of soot clusters formed after the explosion of a hydrocarbon gas mixed with 
oxygen, the cluster morphology was determined by two different methods: structure factor and perimeter 
analysis. Here again we found a hybrid, superaggregate morphology characterized by a fractal dimension 
of D « 1.8 between the monomer size, ca. 50 nm, and 1 pm and D « 2.6 at larger length scales up to 
~10 pm. The superaggregate morphology is a consequence of late stage aggregation in a cluster dense 
regime near a gel point. 
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qa 

Figure 5. — Structure factor, hence expected light scattering pattern, of 3-D simulations of 
aggregation. Late in the run both D = 1.8 and 2.6 morphologies appear. 



Figure 6. — A 1 7 liter chamber 
in which carbon aerosol is 
gelled. 
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Aerosol Gels 


Using the chamber explosions, we discovered a novel one-step method to produce aerogel-like porous 
materials. This method involves the gelation of nanoparticles (primary particles) in the aerosol phase to 
create a material that we have named an aerosol gel. So far we have made both carbon aerosol gels 
(fig. 7) and silica (Si0 2 ) aerosol gels (fig. 8). However, any collection of finely divided primary particles 
with large enough volume fraction can produce an aerosol gel when allowed to aggregate regardless of 
the chemical composition of the parent primary particles. In the case of the carbon aerosol gel the initial 
aerosol is composed of nanometer sized carbon particles produced rapidly by exploding any one of a 
number of hydrocarbons with oxygen in a closed chamber. This method of making an aerosol gel is not a 
wet process and does not require a catalysis. Furthermore, the most important advantage of this method is 
the non-requirement of the supercritical drying step. The carbon aerosol gels have properties quite 
comparable to those of the carbon aerogels. The carbon aerosol gels have high specific area of about 
400 m 2 /g and an extremely low density as low as 2.5 mg/cc. These materials also have high electrical 
conductivity and highly crystalline primary particles. This carbon aerosol gel is significantly different 
than ordinary carbon black and soot formed as byproduct during combustion of any hydrocarbon fuel; it is 
a new material. 

The silica aerosol gel has specific surface areas of ca. 500 m 2 /g. It can be made either hydrophilic or 
hydrophobic depending on the oxidizer used during its synthesis or post gel surface treatment (fig. 9). 



Figure 7. — Carbon aerosol gel. 


Figure 8. — Silica aerosol gel. 



Figure 9. — Carbon coated silica aerosol gel 
with water drop to demonstrate 
superhydrophobicity. 
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Aggregation Kinetics 

Simulations 

Cluster-Cluster aggregation has been studied in our work via large scale Monte Carlo simulations for 
diffusion-limited (DLCA), ballistic-limited (BLCA), and reaction-limited (RLCA) cluster aggregation 
classes in 3-D. We investigated the growth kinetics for systems which evolve from a dilute to a dense 
state, a state we have name “cluster dense,” and then ultimately gel. We studied over three orders of 
magnitude in monomer volume fraction fv. DLCA model results show a power law growth kinetics with 
growth accelerating as the volume fraction occupied by clusters increases with time and the system 
becomes dense. We found that only a weak spatial correlation develops among clusters when dense, and 
cluster-cluster correlation is not sufficient to explain the late time enhancement of the growth kinetics. 
Instead, both exponents, z and X, characterizing the kinetics and size distribution, respectively, show 
universal functionality with cluster volume fraction, independent of the initial monomer volume fraction. 
Remarkably, the relationship between z and X maintains its mean-field nature in 3-D. Results for the 
BLCA model in 3-D is quite similar to the DLCA model while the 3-D RLCA model shows an early time 
exponential growth that crosses over to a power law growth at late times when the system becomes dense. 

We have also studied an off-lattice Monte Carlo simulation of DNA-mediated colloidal assembly. In 
this simulation, aggregation-fragmentation of a binary mixture of DNA-coated colloidal particles is 
studied through a simplified model of base-pair hybridization. Bonding between monomers is modeled as 
a simple temperature sensitive A/B type interaction, where type A and B monomers can only bond to the 
opposite type (no A/A or B/B attachments are allowed). The actual chemistry of base-pair hybridization is 
not included in the model. The morphological structures of the clusters formed as well as the kinetics of 
growth are analyzed in our 2-D simulations. The fractal dimension and kinetic growth exponents for 
clusters formed near the DNA “melting” temperature agree with those seen previously for 2-D diffusion- 
limited cluster aggregation (DLCA) models. The clusters appear more compact, exhibiting signs of local 
order at intermediate temperature values. At higher temperatures, formation of large clusters is not 
favorable under the action of temperature-dependent fragmentation and the system eventually reaches a 
steady state as a collection of small aggregates. The temperature profile for this dissolution of the 
colloidal assembly is sharp, indicating that the selective hybridization process provides a highly sensitive 
measurement tool. At high temperatures, we analyze the steady state behavior of the average cluster size 
in terms of an aggregation-fragmentation model. 

Soot Aerosol Aggregation up to Gelation 

We studied the kinetics, morphology and structure of the quickly aggregating carbon soot particles 
rapidly produced during the explosion of a mixture of a hydrocarbpn gas and oxygen inside a closed 
chamber. We used Small Angle Static Light Scattering (SASLS) technique to do this. Our device allows 
us to detect the scattered light in the angular range O.O5°<0<15°, which corresponds to q values in the 
range 100 cm _1 <q<3xl0“ 4 cm -1 . 

We built a small and thin cylindrical disk combustion chamber with a circular glass window on both 
ends of the cylinder (fig. 10) for the purpose. The internal space of the chamber is 51 wide and 10 mm 
long. The chamber was first evacuated and then filled to different pressures with a hydrocarbon gas- 
oxygen mixture and then ignited using a spark plug. Then the laser light gets scattered by the soot 
aggregates inside the chamber. 

One typical example of the measured scattered light intensity from our small angle SLS device 
plotted versus q is presented in figure 1 1 . The structure factor analysis of the scattered light intensity 
revealed the fractal dimension to be 1.55±0.05. 
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Scattered Light Intensity 



Figure 10. — Small chamber. 
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Figure 1 1 . — Scattered Light Intensity versus q for initial absolute pressure 7 in. of Hg. 
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The Guinier analysis of the scattered light intensity yielded the ensemble average cluster size Rg. 
From Rg versus time we determined the aggregation kernel homogeneity and rate. These are shown in 
figures 12 and 13. Figure 12 shows the homogeneity increases with increasing volume fraction near the 
gel point. The limiting value of X = 0.5 is consistent with our simulations and theory. 




Figure 12. — A plot of the kinetic exponent z and the homogeneity X versus the 
monomer volume fraction f v . The enhanced aggregation kinetics with increasing f v 
was the result of the effect of cluster crowding. 



Figure 13. — A plot K (S 2 , S 2 ) as a function of the monomer volume fraction f v . 
The dashed line shows the theoretical value of K independent of the 
monomer volume fraction. The theoretical K was computed for a 
monodispersed aggregating system in continuum regime with Stokes- 
Einstein diffusion. 
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Aggregation of Charged Aerosols 

We studied the effect of placing charges on a certain fraction of the monomers in an aggregating 
aerosol system. This has relevance to planetary dusts, e.g., Mars, where solar UV photons will charge the 
dust particles via the photoelectric effect, and soot in flames. Maricq has studied the effect of charging on 
the aggregation of carbon soot in ethylene/air laminar diffusion flames. The charges are both positive and 
negative, producing a bipolar population. Due to aggregation, his experimental system eventually reaches 
a point at which the fraction of overall clusters with charges is -33% positive and 33% negative. 

Our study was done through the use of a Monte Carlo algorithm, modified from the standard 
DLCA/RLCA type to include a minimum image type coulomb interaction. The probability of a cluster 
move being accepted during a Monte Carlo time step is based on the Boltzmann factor B = exp(-AE), 
where AE is the difference in coulomb energy difference between the cluster before movement and after. 
The actual probability that a cluster is moved is B/(l+B). Thus, for movements where the energy of the 
final state is significantly higher than the initial, the probability approaches 0. For movements where the 
energy of the final state is significantly lower than the initial, the probability approaches 1 . When there is 
- 0 change in energy, the movement probability is 0.5. Thus there is a deceleration in movements against 
the potential energy gradient and an acceleration in movements in the direction of the energy gradient. 

Simulations were initially performed in the Epstein Regime, appropriate for rarified atmospheres, 
where the motion is still seen as diffusive, but the diffusion constant switches over from the continuum of 
D-Rg” 1 to either D-R g ” 2 or D-A 1 , where A is the cross sectional area of the cluster in the direction of 
cluster movement. Our first simulations were done with the latter, though comparison has been made 
between simulations using all the different diffusion constants. 

For the purpose of our model, monomers are seen as nonconductive, meaning that the charge they 
initially possess is kept throughout the simulation. On aggregation with other monomers/clusters, the 
charges are not spread throughout the cluster. This allows for the possibility of clusters possessing 
nonzero dipole moments. We have also done simulations with conductive monomers, where charges can 
be neutralized on aggregation. Of course, since conductive clusters can contain only one sign of residual 
charge, no cluster dipole moments are possible. For these clusters, the residual charge gets spread among 
the monomers one unit charge at a time in such a way as to minimize the coulomb energy within the 
cluster. 

Specifically, we were interested in seeing the effect of charging on growth kinetics, morphologies, 
and size distributions. As far as kinetics is concerned, the low charge percentage (less than 10% of the 
monomers being charged) used in our simulations does not yield significant differences in the growth of 
clusters. 

In studies of carbon soot in ethylene/air laminar diffusion flames Maricq found that the charge 
distributions broaden as the aggregation proceeds. We observe this in our simulations. However, it does 
not appear that the charges are significantly influencing the aggregation process, as similar results are 
obtained when we turn the charge interactions off and simply label the “charged” monomers. The results 
indicate that charged monomers are simply getting statistically placed in clusters upon aggregation the 
same as any other monomer. This may not be true of systems where the initial population of charged 
particles is higher, and this is something we are continuing to investigate. 

To continue our study of charged aggregation a number of things are being done simultaneously. 

First, as already mentioned, the aggregation of charged systems with a higher initial percentage of charges 
is being investigated both in the case of bipolar and unipolar charging. To better understand the results in 
the Epstein regime, a comparison of the aggregation of uncharged monomers using both possibilities for 
diffusion constant (R g -2 and A -1 ) is being done and a comparison made to the expectations from scaling 
theory. Finally we are developing a Molecular Dynamics program to more correctly model the effect of 
the coulomb force, gravitational settling, etc. 
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Light Scattering 

Mie Scattering 

We continued our previous studies of patterns in Mie scattering. Our tactic to use as an alternative to 
the traditional scattering angle 0 in describing light scattering from a uniform dielectric sphere; that 
alternative is the dimensionless parameter qR, where R is the radius of the sphere, q = 2k sin (0/2) and k 
is the wave number of the incident light. Simple patterns appear in the scattered intensity if qR is used in 
place of 0 that have not been previously recognized (fig. 14). These patterns are characterized by the 
envelopes approximating the scattered intensity distributions and are quantified by the phase shift 
parameter p = 2kR|m-l| where m is the real refractive index of the sphere. We found new patterns in 
these envelopes when the scattered intensity is normalized to the Rayleigh differential cross section. Mie 
scattering is similar to Rayleigh scattering when p<l and follows simple patterns for p>l which evolve 
predictably as a function of p. The curves have a semi-quantitative universality with p. These patterns 
allow us to present a unifying picture of the evolution of Mie scattering for changes in kR and m (fig. 15). 

Multiple Scattering 

We showed that the extent and effect of multiple scattering on angularly resolved light-scattering 
intensity measurements, the optical structure factor, can be quantitatively described by a single parameter, 
the average number of scattering events along the scattering volume. This quantity is easily measured or 
calculated and hence provides a useful experimental indicator of multiple scattering, which is a hindrance 
to accurate structure factor measurements. 



Figure 14. — Scattered light intensity normalized by the Rayleigh cross section for 
spheres of arbitrary radius R and refractive index m versus qR. Note quasi 
universality on the shift parameter p = 2kR|m-1 1. 
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